Unified approach to QED in arbitrary linear media 
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We give a unified approach to macroscopic QED in arbitrary linearly responding media, based 
on the quite general, nonlocal form of the conductivity tensor as it can be introduced within the 
framework of linear response theory, and appropriately chosen sets of bosonic variables. The for- 
malism generalizes the quantization schemes that have been developed previously for diverse classes 
of linear media. In particular, it turns out that the scheme developed for locally responding linear 
magnetodielectric media can be recovered from the general scheme as a limiting case for weakly spa- 
tially dispersive media. With regard to practical applications, we furthermore address the dielectric 
approximation for the conductivity tensor and the surface impedance method for the calculation of 
the Green tensor of the macroscopic Maxwell equations, the two central quantities of the theory. 

PACS numbers: 42.50.Nn, 03.70.+k, 12.20.-m 



I. INTRODUCTION 

In both classical and quantum electrodynamics, it is 
often advisable to divide, at least notionally, the matter 
that interacts with the electromagnetic field into a part 
that plays the role of a passive background and a re- 
mainder, active part that needs to be considered in more 
detail. By means of suitable coarse-graining and aver- 
aging procedures, this leads to the well-known frame- 
work of Maxwell's phenomcnological equations, where 
the background — the medium — is treated as a continuum 
and, quite frequently, by the methods of linear response 
theory. From this perspective, the characterization of the 
medium is reduced to the prescription of suitable consti- 
tutive relations, i.e., appropriate response functions or 
susceptibilities. 

Depending on the specific kinds of media under con- 
sideration, it is under many circumstances sufficiently ac- 
curate to work with spatially local response functions, 
taking into account only (temporal) dispersion and ab- 
sorption in accordance with causality. For conducting 
and semiconducting media (not to mention plasmas) as 
well as superconducting materials, however, the spatially 
local description can be inadequate due to the exis- 
tence of almost freely movable charge carriers (conduc- 
tion electrons, excitons, Cooper pairs) in such media. 
Hence, if one is not willing to restrict one's attention to 
a crude spatial resolution and/or specific frequency win- 
dows, spatial dispersion, i.e., the spatially non-local char- 
acter of the medium response, generally cannot be disre- 
garded for such media. Electrodynamics problems with 
the inclusion of spatial dispersion have been considered 
by various authors in different ways, both on the classi- 
cal and quantum levels; for classical approaches, see, e.g., 
Refs. [1—7], for quantum ones see, e.g., Refs. [8, 9]. 

A scheme that takes spatial dispersion into account 
along with dispersion and absorption in sufficiently gen- 
eral terms can also be regarded as an important step to- 
wards a satisfactory (quantum) electrodynamics of mov- 



ing media, which is very much lacking at present. The 
reason is that a medium, even if it can be assumed to re- 
spond spatially locally when it is at rest, will in general 
appear as responding non-locally when it is in motion. 
Given that the polarization of typical Drude-Lorentz- 
type dielectrics responds to the electric field with a char- 
acteristic memory time of the order of 10~ 9 ■ ■ • 10~ 7 s 
[10], already moderate (i.e., non-relativistic) velocities 
may lead to the appearance of noticeable spatial non- 
localities. For example, sonoluminescence experiments 
show that the collapse of a bubble with a typical initial 
radius of 10 • • • 50/im to a final radius of around 1/itm oc- 
curs on a time scale similar to the characteristic memory 
time of the response of the surrounding fluid [11]. 

The study of the quantized electromagnetic field in 
spatially non-locally responding media and the prospect 
of elaborating a quantum theory of light in moving media 
will also open up new ways of investigating quantum ef- 
fects related to the recently proposed 'optical black hole' 
[12, 13]. So far, the theory has concentrated on purely 
geometrical optics with some progress being made to- 
wards a (scalar) wave-optical description, but a consis- 
tent linear-response approach is still lacking. 

The paper is organized as follows. In Sec. II we intro- 
duce the basic concepts of field quantization in arbitrary 
linearly responding media, with special emphasis on spa- 
tially dispersive media. This serves as the basis for a 
detailed study of possible choices of appropriate dynam- 
ical variables in Sec. III. We then proceed to show in 
Sec. IV how previously introduced quantization schemes 
for diverse classes of media can be obtained as special 
cases from the general quantization scheme developed in 
Sec. II. In addition to the general formalism, some knowl- 
edge of the structure of the Green tensor for spatially 
dispersive media is needed when performing explicit cal- 
culations. This problem is addressed in Sec. V, where 
it is described (in general and by an example) how the 
surface impedance method may be applied in this con- 
text, on the basis of the dielectric approximation. Some 
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concluding remarks are given in Sec. VI. 



II. QUANTIZATION SCHEME 

The effect of any linear, dispersing and absorbing 
medium on the electromagnetic field can be described, 
within the framework of linear response theory, by the 
relation 

j(r,w) = |dVQ(r,r', W ).E(r', W )+j N (r, W ), (1) 

where j(r, u>) and E(r, u>), respectively, are the (linearly 
responding) current density and the electric field in the 

frequency domain, Q(r, r',ui) is the complex conductiv- 
ity tensor in the frequency domain [7, 14], and j N (r, ui) is 
a Langevin noise source. According to the Onsager reci- 
procity theorem [7, 14] , the conductivity tensor should be 
reciprocal, Qij(r, r', ui) = Qji(r' , r, ui). Except for a trans- 
lationally invariant (bulk) medium, the spatial arguments 

r and r' of Q(r, r', u) must be kept as two separate vari- 
ables in general. We assume that, for chosen u, Q(r, r',ui) 
is the integral kernel of a reasonably well-behaved (in- 
tegral) operator acting on vector functions in position 

space. In particular, we assume that Q(r,r',ui) tends 
(sufficiently rapidly) to zero for |r — r'\ — > oo and has no 
strong (i.e., non-integrable) singularities (specifically, for 
r' — > r). To allow for the spatially non-dispersive limit, 
(^-functions and their derivatives must be permitted so 

that Q(t,t',u>) may become a (quasi-)local integral ker- 
nel. In the remainder of the paper, we will use the su- 
perscripts T and + to indicate transposition and Hcrmi- 
tian conjugation with respect to tensor indices. Since the 
spatial arguments are not switched by these operations, 

the operator associated with an integral kernel A(r, r') is 

Hermitian if A(r, r') = A + (r' , r). In particular, an opera- 
tor associated with a real kernel is Hermitian if it has the 
reciprocity property A(r, r') = A J (r' , r). The decompo- 
sition Q(r, r',ui) = Re Q(r, r',ui) + i Im Q(r, r', ui) of the 
conductivity tensor is therefore identical with the decom- 
position of the associated operator into a Hermitian and 
an anti-Hermitian part, 

a(r,r',u) = ReQ(r, r',ui) 

= ±[Q(r,r>) + Q + (r',r,^)], (2) 

r (r, r',ui) = Im Q(r, r', u) 

= ±[Q(r,r',u>)-Q+(r , ,r,u>)]. (3) 

Since a (r, r', u>) is associated with the dissipation of elec- 
tromagnetic energy (see, e.g., Refs. [7, 14]), the operator 
associated with the integral kernel a(r,r',uj) is, for real 



ui, a positive definite operator in the case of absorbing 
media considered throughout this paper. 

The conductivity tensor Q(r, r',u>) is the temporal 

Fourier transform of a response function Q(r, r', t) in the 
time domain, 



Q(r,r» = J dte^Q(r,r',t), (4) 
which satisfies causality conditions of the type 



Q(r,r',t)=0 if t 



cos r\ |r 



|r-r'|/c<0 (5) 



for chosen r and r' and arbitrary directional cosines cos r] 
(0 < cos?] < 1). In particular, for cos?7 = 0, one finds 
from arguments [14-16] similar to those for the case 
of spatially locally responding media that, for chosen r 

and r', Q(r, r',u>) is analytic in the upper complex u> 
half-plane, fulfills Kramcrs-Kronig (Hilbcrt transform) 
relations, and satisfies the Schwarz reflection principle 

Q*(r,r',u)) — Q(r, r', — u>*). Other values of cos i] could 
obviously provide more stringent (spatio-temporal) con- 
ditions (see also Ref. [17]), which are, however, not re- 
quired here. 

Let us identify the current density that enters the 
macroscopic Maxwell equations in the frequency domain 
with j(r, ui) as specified in Eq. (1). In this case, the 
medium- assisted electric field in the frequency domain 
satisfies the integro-differential equation 

V x V x E(r, u) tt E(r, w) 

c l 

- ifi u> / d 3 / Q(r, r', w) ■ E(r', u>) = i/j, u>j (r, w), (6) 



whose unique solution is 



E(r, u>) = i» u> / dV G(r, r', w) • j fr', w), (7) 



with G(r, r',ui) being the (retarded) Green tensor. It 
satisfies Eq. (6) with the (tensorial) 5-function source, 



V x V x G(r, s,ui) 2" G{r, s, uj) 



— IjlQU) 



Q(r,r',w) -G(r',s,uj) = /<5(r-s), (8) 



together with the boundary condition at infinity, and 
has all the attributes of a (Fourier transformed) causal 

response function just as Q(r,r',u>) has them. In par- 
ticular, it is analytic in the upper u> half-plane and the 

Schwarz reflection principle G*(r, r',ui) = G(r, r', — u>*) 
is valid. Its basic properties in position space are simi- 
lar to the ones known from the spatially local theory, in 

particular, since Q(r, r',ui) is reciprocal, so is G(r,r',w), 
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G(r, r', ui) = G T (r', r, ui), and, for real uj, the generalized 
integral relation 

pouJ J d 3 s J d 3 s'G(r,s,w) • a(s, s',uj) • G*(s',r',w) 

= ImG(r,r» (9) 

holds (App. A). 

To quantize the theory, the Langevin noise source 
3 N (r,oj) is regarded as an operator |j N (r,w) i-> j N (r,w)] 
with the commutation relation 



[j (r.wJjJ^rV)] = —5{ui-uj')a{v,v l ,uj). (10) 

The (now operator-valued) equation (7) relates the elec- 
tric field operator 

/>OC 

E(r) = / dwE(r,w) + H.c, (11) 
Jo 

and thus all the electromagnetic field operators, to 
j N (r, uj) and jjj(r',a/), which may be regarded as the 
dynamical variables of the overall system consisting of 
the electromagnetic field and the linear medium (incor- 
porating the reservoir degrees of freedom responsible for 
absorption). It should be mentioned that, by means of 
the correspondence 



e uj 



Q(r,r',u>) <-» x(r,r',w), 



(12) 



where x( r > r ': w ) is the (nonlocal) dielectric susceptibil- 
ity tensor, the basic commutation relation (10) becomes 
equivalent to the commutation relation derived from a 
microscopic, linear two-band model of dielectric material 
[18], which has been used to study the quantized electro- 
magnetic field in spatially dispersive dielectrics [8, 9]. 

In order to complete the quantization scheme, a Hamil- 
tonian H needs to be introduced [as a functional of 
j N (r,w) and jjj(r', uj')] so as to generate 'free' time evo- 
lution according to 

g N (r,w),fl] = fiwj N (r,w) ) (13) 
which constrains the Hamiltonian to the form 

H = irj°° J d 3 r J dV j^r, w) • p (r, r', uj) • j N (r', uj), 

(14) 

to within irrelevant c-number contributions. A glance 
at Eqs. (10) and (13) now shows that p(r,r',ui) is the 
integral kernel of the inverse operator of the operator 
associated with cr(r, r' ,ui) (which exists). The validity 
of the quantization scheme is confirmed by checking that 
the well-known equal-time commutation relations for the 
electromagnetic field operators hold, which can be done, 
in analogy to the spatially local theory (cf. Refs. [19- 
21]), by properly taking into account the properties of 

G(r, r',ui) and Q(r, t' ,ui) [in particular, Eq. (9)]. 



The Hamiltonian (14) may clearly be brought to the 
diagonal form 

H = Jd 3 rJ^ dujfiuj^(r,uj)-f(r,uj) (15) 

known from the spatially local theory, where f (r, uj) is a 
bosonic field, 

[f(r,a,),ft(r>')] = 6{w - u/)I6{r -t>), (16) 

by performing a linear transformation of the variables, 
which we shall assume to be invertible. Writing 



j N (r,c)=^— j jdV^(r,r».f(r', W ), (17) 

the diagonalization is achieved and Eqs. (16) and (10) 
are rendered equivalent if we choose the integral kernel 

K(r, r',u) such that, for real ui, the integral equation 



d 3 sif(r,s,w) • K + (r',s,oj) = a{r,r',u) (18) 



holds, which is guaranteed to possess solutions (sec 
Sec. Ill) since <r(r, r',uj) is the integral kernel of a posi- 
tive definite operator. 

So far we have considered the 'free' medium-assisted 
electromagnetic field. Its interaction with additional 
(e.g., atomic) systems can be included in the theory on 
the basis of the well-known minimal or multi-polar cou- 
pling schemes in the usual way (see, e.g., Ref. [22]). 



III. NATURAL VARIABLES AND PROJECTIVE 
VARIABLES 

Let us now turn to the problem of constructing the 

integral kernel K(r,r',u)) in Eq. (18). For this purpose, 
we consider the eigenvalue problem 

J dV a(r, r', uj) ■ F(a, r',u) = a(a, w)F(a, r, u) (19) 

which, under appropriate regularity assumptions on the 

conductivity tensor Q(r, r',u), such as those listed be- 
low Eqs. (1), is well-defined. In particular, it features a 
real (positive) spectrum and a complete set of orthogonal 
eigensolutions, which we may take to be (J-)normalized. 
Note that the real u> plays the role of a parameter here, 
and a stands for the collection of (discrete and/or contin- 
uous) indices needed to label the eigenfunctions. Adopt- 
ing a continuum notation, we may write 

J daF(a,r,w)F*(a,r»= ?S(r-r'), (20) 
J d 3 rF*(a,r,w).F(a',r,u)=i(a-a'), (21) 
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and the diagonal expansion of a (r, r', u) reads 

<r(r, r', w) = j da a (a, u)F(a, r, uj)F* (a, r', w), (22) 

which resembles the expansion of the dielectric suscep- 
tibility in the microscopic theory [18] mentioned above. 
Substituting Eq. (22) into Eq. (18), we may construct an 

integral kernel K(r, r', ui) in the form of 
K(r,r',w) = ( dacri(a,uj)F(a,r,Lu)F*( a ,r , ,uj), (23) 



where we choose ct 1 / 2 '(a, u) > so that the operator asso- 
ciated with K(r,r' ,u>) is the positive, Hcrmitian square- 
root of the operator associated with cr(r, r',w). Obvi- 
ously, this solution to Eq. (18) is not unique, since any 
other kernel of the form 

K'(r,r',w)= J d 3 sK(r,s,w)-V(s,r',w) (24) 

with V(r,s,u>) satisfying 

d 3 sV(r,s,u) •V r+ (r / ,s,w) = Ts(r - r') (25) 



also obeys Eq. (18). As we are interested in invert- 
ible transformations (17), the operator corresponding to 

V(r, s, u) should be invertible as well, so that we can 
replace Eq. (25) with the stronger unitarity condition 



/ 



d 3 sV + (s,r,uj)-V(s,r',u}) 

= I d 3 s V(r, s, u>) ■ V + (r' } s, w) = Ts(r - r'). (26) 



/■ 

Without loss of generality (see App. B), we can base our 
further calculations on Eq. (23). 

Inserting Eq. (23) into Eq. (17), we find that 

J N ( r ' w )= J daa^(a,u;)F(a,r,uj)g(a,u}) 7 

(27) 

where we have introduced the new variables 

g(a,w) = J d 3 rF*(a,r,co) ■ f( T ,w) , (28) 

referred to as the natural variables in the following. 
Needless to say that they are again of bosonic type, 

[g(a,w),g\a',w')] = 5 (a - a/)5(u> - J). (29) 

Since the transformation (28) does not mix different 
lo components, the Hamiltonian (15) is still diagonal 
when expressed in terms of the natural variables, 

H = J da J duj fuv g^ (a,uj)g(a,uj), (30) 



as can be easily seen by inverting Eq. (28), 

f(r,w) = / daF(a,r,Lu)g(a,Lu), (31) 



and combining with Eq. (15), on recalling Eq. (21). 

Let us organize the set of eigenfunctions F(a,r, u) 
into (a discrete number of) subsets labeled by A (A — 
1,2,..., A). With the notation a i — > (A,/3),Eq. (31) then 
reads 

f(r,w) = ^f A (r,u;), (32) 

A 

where 

f*(r,w) = J d(3F x {(3,v,u)g x {p,uj). (33) 
The operators associated with the integral kernels 

P A (r,r» = / d/3F A (/?,r,u;)F A (/?,r>) (34) 



form a complete set of orthogonal projectors. Obvi- 
ously, these projectors and the operators associated with 

a(r, t',uj) and K(r, r',w) as given by Eq. (23) are com- 
muting quantities. It is not difficult to see that the vari- 
ables 

f\(r,u) = J dVP A (r,r', W ).f(r', W ) 

d/?F A (/3,r,w)&(/?,w), (35) 



referred to as projective variables in the following, obey 
the non-bosonic commutation relation 

[f A (r,w),ft,(r>')] =5 xx ,S(w-u')P x (r,r',u,), (36) 

and the Hamiltonian (15) expressed in terms of the pro- 
jective variables reads as 

/roc 
d 3 r J dwfcjft(r,w)-f A (r,w). (37) 

From Eqs.(36) and (37) it then follows that 

[f x ( T ,u),H] = hwJd 3 r'P x {v,r',u)-f x {r',u) 

= hwf x (r,u>). (38) 
Inserting Eq. (32) in Eq. (17), we obtain 

A 

where the J NA ( r , w) are given by 
1(t,w)=(—Y fdVf A (r,r».f A (r», (40) 



7T / 
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with 

K x (r,r',oj) = J d 3 S P x (r,s,uj)-K(s,r',oj) 

= J d 3 S K(r,s,Lu)-P x (s,r',u). (41) 

Recalling Eq. (36), we can easily see that 

[l NA (r,^)JNA'( r 'V)] = ^6x\>6(w-w')cTx(r,r',w), 

(42) 

where <j\(r,r',uj) is defined according to Eq. (41) with 

a in place of K. Summation of Eq. (42) over A and 
A' leads back to Eq. (10), so that the two equations are 
equivalent. 

At this stage, we observe that there is the option to 
base the quantization scheme directly on Eqs. (37), (39), 
and (40), regarding the variables f\(r, w) and f|(r, oj) as 
the basic dynamical variables of the theory and assigning 
to them bosonic commutation relations 

[f A (r, W ),f],(r>')] =5xx'S(u-uj , )Ts(T-r') (43) 

in place of Eq. (36). Note that, in so doing, back refer- 
ence from the variables f\(r,u>) to the original variables 
f(r,a;) is not possible anymore. As can be seen from 
Eqs. (40) and (41), Eq. (42) is satisfied also when the 
f\(r, oj) and f|(r, oj) are considered as bosonic variables, 
from which it follows [via Eq. (39)] that Eq. (10) also 
still holds and, as before, this implies that the correct 
electromagnetic-field commutation relations hold. The 
second line of Eq. (38) remains of course also true so 
that the correct time evolution is ensured as well. 

Since the state space attributed to the bosonic vari- 
ables f\(r, oj) and f|(r,w) is, in general, different from 
the state space attributed to the original variables f (r, oj) 
and F(r, oj) [or, equivalcntly, attributed to g x ((3,0j) and 
g\{(3,uj)], the allowable states must be restricted, by rul- 
ing out certain coherent superpositions of states in the 
sense of a super-selection rule. In App. C, we show that 
the condition imposed on the states may be described by 
means of a set of projectors P x such that the allowable 
states can be characterized by 

P A |V> = |V)VA, (44) 

where the action of the projectors P\ in state space is 
closely related to the action of the projectors associated 
with the kernels (34) in position space. As a result, if 
the total Hamiltonian H to t composed of the Hamiltonian 
(37) and possible interaction terms (in the case where 
additional, active sources are present) commutes with all 
of the projectors Px, 

[Px,Htot] =0VA, (45) 

then allowable states remain allowable in the course of 
time, and the option of treating the f\(r, w) and f|(r, oj) 



as bosonic variables can be safely exercised. Clearly, 
all the observables of interest should then also commute 
with the Px so that no transition matrix elements be- 
tween states belonging to different subspaces, i.e., be- 
tween spaces attributed to different A values, can ever 
come into play. 

One can also consider decompositions of j N (r, oj), 
where in place of the j NA (r,c<j) introduced above other 
quantities J NA (r, uj) subject to the condition 

^J NA (r,c) =^j NA (r,c) (46) 

A A 

are introduced, whose commutation relations may be 
quite different from those of the j NA (r, w). Obviously, 
the total noise current density j (r, w) as defined by 
Eq. (39) and the commutation relation (10) arc not 
changed by such a transformation, briefly referred to as 
gauge transformation in the following. Moreover, since, 
with regard to Eq. (10), only the sum of the commu- 
tators [J NA (r, uj), (r', uj')] over all A and A' is rele- 
vant, every chosen set of (algebraically consistent) com- 
mutators [J NA (r, oj), jJf A /(r', oj')] which leads to Eq. (10) 
yields, in principle, a consistent quantization scheme in 
its own right. A 'substructure below' Eq. (10) can hence 
be introduced with some arbitrariness, but since the var- 
ious available alternatives are not necessarily equivalent 
to each other, a specific one should not be favored in 
the absence of good (physical) motivation. In contrast, if 
the observables of interest — including the Hamiltonian — 
can be viewed as functionals of J N (r, oj) [rather than of 

the individual J NA (r, w)], Eqs. (10) and (14) can be re- 
garded, in view of the fluctuation-dissipation theorem(s) 
(see, e.g., Ref. [14]), as being unique, and hence, as in- 
variable fundament of the theory. 

From the above, it may be reasonable to widen the 
notion of projective variables as follows. If, for a chosen 
(physically motivated) decomposition of the noise current 
density, it is possible to (linearly) relate the J NA (r, w) in 
Eq. (46) to (new) variables f\(r,cj) such that, upon con- 
sidering the latter as bosonic variables, the validity of the 
basic equations (10) and (13) is ensured, then the spe- 
cific quantization scheme so obtained may be regarded 
as arising from the general quantization scheme by ex- 
cluding certain types of (superposition) states from state 
space, and restricting the dynamics (as well as the allow- 
able observables) accordingly. The f \ (r, oj) may then be 
seen as projective variables in a wider sense. 



IV. DIFFERENT CLASSES OF MEDIA 

We proceed to show that rather different classes of me- 
dia (usually studied separately) fit into the general quan- 
tization scheme developed in Sec. II. The main task to be 
performed is solving the eigenvalue problem (19), which 
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requires knowledge of a (r, r', to) for the specific medium 
under consideration. In two limiting cases, the exact so- 
lution to Eq. (19) can be given straightforwardly, namely, 
in the case of an inhomogeneous medium without spatial 
dispersion and in the case of a homogeneous medium that 
shows spatial dispersion. Let us, therefore, first examine 
these two cases in detail before considering more general 
situations. 

A. Spatially non-dispersive inhomogeneous media 

The complete neglect of spatial dispersion means to 

regard the medium response, i.e., Q(r, r',uj), as being 
strictly local. If this is assumed, we have 

a(v,v',uj) = a(v,uj)6(r-v'), (47) 

where a (r, uj) can be written in diagonal form as 

3 

^(r,^) = ^<^i{'r^)e l {v,uj)e*{v,uj), (48) 

i=l 

with ej(r, uj) (i = 1,2,3) being orthonormal unit vec- 
tors. Hence, the eigenvalues a(a,u>) and eigenfunctions 
F(a, r, uj) of the operator associated with a(r, r',oj) read 
[a i ► (i,s)] ai(s,u)) and 

Fj(s,r, uj) = ei(s,uj)5(s - r), (49) 

respectively. Equation (23) then becomes 

K(r,r',uj) = K(r,uj)6(r-r'), (50) 

where 

^ 3 

K{y,uj) = ^cr l 1/2 (r,w)e,(r,w)e*(r,w), (51) 
»=i 

and Eq. (17) takes the form 

j N (r, W )=(^) 2 f(r, W )-f(r, W ), (52) 

which just yields the well-known quantization scheme for 
a locally responding, possibly anisotropic dielectric ma- 
terial [19, 20], upon identifying a(r, uj) = Sou)lm x(r, uj), 
with x(r, uj) being the (local) dielectric susceptibility 
tensor [cf. Eq. (12)]. The natural variables §i(r, uj) are 
here simply the components of f(r, uj) along the princi- 
pal axes of the medium, which may in general vary with 
position and frequency, 

gi(r, uj) = e*(r, uj) ■ f (r, uj), (53) 

3 

f ( r , w ) = e *( r ' w )&( r ' w )- ( 54 ) 
»=i 



Identifying the index A introduced in Eq. (32) with i 
and assuming that <7j(r, uj) ^ (r, uj) for i^i', one can 
define, according to Eq. (34), the three projection kernels 

Pi(r,r',uj) =e i (r,w)e?(r / ,a;)5(r-r / ), (55) 

which, according to Eq. (35), give rise to three sets of 
projective variables, 

f»(r, uj) = e,(r, uj)gi(r, uj). (56) 

As long as the projective variables are not coupled to each 
other — which is obviously the case for the 'free' system 
governed by the Hamiltonian (38) — they can be regarded 
as being of bosonic type. In this case, instead of using 
the original set of bosonic variables f(r,uj) and r(r,uj), 
one can use three sets of bosonic variables fj(r, uj) and 
{■(r,uj) associated with the three principal axes of the 
dielectric medium at each space point. 

If two of the three eigenvalues <7j(r, uj) coincide (uni- 
axial medium), the two corresponding projection kernels 

Pi(r,r',w) should be combined into one projector (pro- 
jecting on the plane perpendicular to the distinguished 
axis of the medium), thereby reducing the number of sets 
of projective variables to two. Clearly, if the three eigen- 
values <7j (r,u) all coincide (isotropic medium), the three 

projection kernels Pj(r, r', uj) should be combined to give 

the unit kernel 1 5(r — r'), corresponding to the use of 
the original variables. 

B. Spatially dispersive homogeneous media 

In the limiting case of an (infinitely extended) homo- 
geneous medium, Q(r,r',uj) is translationally invariant, 
i.e., it is a function of the difference r— r', and so is then 
a(r,r',uj). We may therefore represent it as the spatial 
Fourier transform 

a(r, v',uj) = ^ J d 3 k a (k, uj)^'^ , (57) 
where 

3 

<r(k,w) = <7j(k, uj) ej(k, u;)e*(k, uj), (58) 

i=l 

with ej(k, uj) (i = 1,2,3) being orthogonal unit vectors. 
Consequently, the eigenvalues <r(a, uj) and eigenfunctions 
F(a,r,uj) of the operator associated with <r(r, r',uj) are 
[a i ► (i,k)] <Ti(k,uj) and 

Fi(k, r, uj) - (2^)- 3 /V k - r e 4 (k, uj), (59) 
respectively, and Eq. (23) reads 

K{v,v',uj) - ^ J d 3 kK (k, W )e ik '( r - r '), (60) 
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where 



K(k, W ) = ^<7, 1/2 (k, W )e i (k,a;)e*(k, W ). (61) 



Combination of Eqs. (17), (60), and (61) then yields 



J 1 v( r > w ) = ~ 



Eq. (18), since any other kernel K'(r,r',uj) of the form 
(24) [together with Eq. (26)] is also an allowed kernel. To 
illustrate this for the isotropic medium under study, we 
first note that Eq. (65) may be equivalently rewritten as 



a (k, uj) = <7|| (k,uj)I - k x j(k, w)I x k, (69) 



where 



V 7T J (2tt) 3 / 2 



7(fc, LO) = [(7±(fc, Uj) — (Til (fc, U))]/k 2 . 



(70) 



r 

xJ2 d 3 fce lk - r a 2 1/2 (k,^)e 4 (k, W )5 4 (k, W ), (62) 
i=i 

where the natural variables gi(k, uS) are related to the 
spatial Fourier components of f (r, uj) as 



g^) = -^^ J d 3 re-^e*(k,u)-f(r,u). 



(63) 



On the basis of the three unit vectors ej(k, w), three 
(different) projection kernels can be introduced, 

P i (r,r» = 7 ^ / d 3 fce 4 (k,^)e*(k, W )e tk -( r - r '), 



(27T) 3 . 

(64) 

provided that <7,(k, w) 7^ ay (k, uj) for i'. 

Let us consider, in particular, isotropic media that in 
addition do not feature optical activity in more detail. In 
this case, the diagonal form of the tensor a (k, uj) reads 
(see Ref. [7]) 

«-> kk kk\ 

a(k,uj) = a\\(k,uj)— + aj_(k,uj) (/-_), (65) 



Since (for real uj) a\\{k,uj) and oj_(fc,cj) are both real 
and positive [in accordance with the requirement that 
a (r, t',w) be the integral kernel of a positive definite op- 
erator], 7(fc,w) is real but its sign is not determined by 
this requirement. However, if j(k,u>) is required here and 
below to be positive throughout, then 

K'(k,u) =a\ /2 {k,uj)7 ± 7 1 / 2 (/c,w)k x F (71) 

obeys the equation 

K'{k,ui) ■ K'+{k,uj) = a{k,uj). (72) 

Moreover, it can be shown that K'(k, ui) can be repre- 
sented in the form 



K'(k,uj) = K(k,uj) -V(k,uj), 



(73) 



with 



nk,u)^ + o\'\k,u)a^\k,u)(7 

± 7 1 / 2 (fc, Lu)a~ 1/2 (k, Lu)k x J (74) 



i.e., o-i(fc,w)=CT||(fc,a;) and a 2 (fc, w) = CT 3 (fc, u>) =a±(k, u>) [V' 1 (k,uj) = V + (k,ui)]. Hence, K'(k,ui) also yields, ac- 



^(7||(fc,cj), which implies that K(k,uj), Eq. (61), takes 
the form 

K(k, uj) = a\ /2 {k, u)^+ ai /2 (fc, u )(f-^j. (66) 

Thus, the well-known longitudinal and transverse tenso- 

*— * <— » 

rial (5-functions A||(r — r') and Aj_(r — r'), respectively, 
can be taken as projection kernels, 



cording to Eq. (60), a valid integral kernel K'(r, r',uj), 
£'(r,r» = -^3 J d 3 kK'(k,u;)e ik < r - r '\ (75) 

which is related to the integral kernel K(r,r',uj) accord- 
ing to Eq. (24), where 



y(r,r>) = ^3 J d 3 fcy(k,o;) e ik -( r - r '), 



(76) 



P u±) (r,r',uj) = A| K±) (r-r'), 



(67) 



which may be used to introduce, according to Eq. (35), 
the projective variables 



f||(±)(r,w) 



(r-s) -f(s,w). (68) 



Unitarily equivalent formulation 

As already pointed out in Sec. II, the integral ker- 
nel K(r,r' ,w) in Eq. (17) is not uniquely determined by 



with the associated operator being unitary. We thus 
see that the two formulations of the theory based on 

K(r, r',w) and ^'(r, r', w), respectively, are unitarily 

equivalent. Note that K'(k,uj) ^ K' + (k,uj), so that the 

operator associated with the integral kernel K'(r,r',uj) 
is non-Hermitian (as to be expected, see App. B). Since 

the operators associated with K'(r,r' ,u>) [as well as 

V"(r, r',uj)} and P||(j_)(r, r', ui) commute, the same pro- 
jectors may be employed in the two formulations of 
the theory to introduce projective variables according to 
Eq. (68). 
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2. Local limit: magnetodielectric media 

Now let us suppose that <j\\{k,uj) and -f(k,uj) in 
Eq. (69) are sufficiently slowly varying functions of 
k, with well-defined and unique long-wavelength lim- 
its linife^o (?\\{k,uj) = <t\\{lu) > and lim^o l{k, w) 
= 7(0;) > 0, so that they may be approximated by these 
limits under the integral in Eq. (57) to obtain 

a(r,r',uj) = ai\(u>)IS(r - r') 

-7(w)Vx[H(r-r')]xV'. (77) 

It should be pointed out that in the limiting case given 
by Eq. (77) the positive definiteness of the operator as- 
sociated with tr(r, r',u>) already implies that j(ui) must 
be positive, 7(0;) > 0; in the general case as given by 
Eq. (57) together with Eqs. (69) and (70), the positive 
definiteness of a (r, r',w) does not automatically restrict 
j(k,u>) to positive values. 

In order to see to what type of medium this cr(r, r' ,ui) 
corresponds, we have to find from Eq. (77) the full con- 
ductivity tensor Q(r,r',u), which is uniquely possible 
since Eqs. (2) and (3) are Hilbert transforms of each other 
(cf. Sec. II). The full conductivity tensor corresponding 
to Eq. (77) is thus of the form 

Q(r,r',w) = Q( 1 )( W )^(r-r') 

-Q (2) (co)V x [/J(r-r')] x V', (78) 

where Q^(u>) and Q^ 2 \lo) are (Fourier-transformed) re- 
sponse functions, both of which are determined by their 
respective real parts (T\\{uj) and 7(0;). Inserting Eq. (78) 
into Eq. (1) and comparing with 

j(r,w) = -ie a uj [e(ui) - l]E(r,w) 

+ k V x {[l-K(w)]B(r,o;)}+j N (r,a;) (79) 

[B(r,w) = (iw) _1 V x E(r, w)], which is the well-known 
description of a locally responding (homogeneous) mag- 
netodielectric medium, we can make the identifications 

Q^(u) = -is u)[e(u)-l} (80) 

and 

Q {2 \lo) = -ino[l - k(u>)]/u>, (81) 

where s(lo) is the permittivity and n(to) = k (w) the 
(paramagnetic) permeability of the medium ((j,q = k,q ). 
For real u, we thus obtain 

<7\\(oj) = £owIme(w) (82) 

and 

7(0;) = — Kolm k(u>)/u>. (83) 



Note that, because of 7(0;) > 0, from Eq. (83) it follows 
that Im/t(o;) < for u> > 0, from which it can be shown 
that >0) > 1. 

At first glance, one might believe (erroneously) that 
not only paramagnetic [fj,(oj — > 0) > 1] but also diamag- 
netic [/i(o-> — > 0) < 1] features of a medium (or the com- 
bined effect of both) can be consistently described by 
means of the magnetic permeability /i(uj) which is in- 
cluded, as seen above, in the basic linear-response con- 
stitutive relation (1). However, since diamagnetism is ba- 
sically a nonlinear effect (as the underlying microscopic 
Hamiltonian is quadratic in the magnetic induction field), 
it is beyond the scope of linear response theory. If it is 
desired to include diamagnetic media in the framework of 
linear electrodynamics nevertheless, one can regard the 
magnetic field on which the diamagnetic susceptibility 
depends as being (the mean value of) an externally con- 
trolled field independent of the dynamical variables. Note 
that the Onsager reciprocity theorem needs to be stated 
in its generalized form in this case, see Refs. [7, 14]. For 
a more satisfactory account of diamagnetic media, one 
should, however, resort to a non-linear response formal- 
ism, or to a more microscopic theory. 

An obvious solution to Eq. (18) with <r(r, r',w) given 
by Eq. (77) is provided by 

K'(r,r',(j) = o- | ] /2 (a))F(5(r-r')Tj7 1/2 (^)Vx Ts(r-r'), 

(84) 

which corresponds to the kernel (75) [together with 
Eq. (71)] when, for an isotropic medium, a\\{k,uj) and 
j(k,u) are approximated by (T\\(u)) and respec- 
tively. The kernel (75) [together with Eq. (71)] fits well 
here since it depends in a particularly simple way on 
those quantities that we have assumed to approach well- 
defined limits in the derivation that led to Eq. (77), 
a property which can be attributed to the responsi- 
ble transformation (76) [together with (74)]. In con- 
trast, the kernel obtained directly from Eq. (60) [together 
with Eq. (66)], by first eliminating a±(k,u) by means 
of Eq. (70) and then approximating au(k,u>) 1— > <7||(u;) 
and j(k,u)) 1— >7(w), does not provide an alternative to 
Eq. (84), as it does not lead to Eq. (77) when inserted 
in Eq. (18); it corresponds to a different medium. Corre- 
spondingly, the kernel V(r, r',w) obtained by expressing 
Eq. (76) [with Eq. (74)] in terms of a\\(k,u>) and j(k,uj) 
and then approximating them by aii (u) and j(cj), respec- 
tively, fails to be unitary. 

Substituting for K(r, r',w) in Eq. (17) K'(r,r',u) as 
given by Eq. (84) [together with Eqs. (82) and (83)], we 
may explicitly express the noise current density in terms 
of the bosonic dynamical variables to obtain 

1 

=F i (—) 2 V x [y/-hn k{uj) f (r, w)] . (85) 
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Since the operators associated with the projection ker- 
nels (67) commute with the operators associated with 
Eqs. (84) and (77), one may introduce the projective 
variables fi|(j_)(r, w) defined by Eq. (68), which corre- 
sponds to a decomposition of j (r,w) into longitudinal 
and transverse parts, 



j N (r,^) = j N|| (r,w) +j N± (r,w), (86) 



where 



he 



j N|| (r, W )=(^) V^ImeMfiiM, (87) 



he 



\/u 2 Ime(w) fj_(r, ui) 



Tl {^r) Vx [V-Im«(w)f ± (r,w)]. (88) 

Making use of Eq. (35) and identifying therein the projec- 
tion kernels P\(r, r', ui) with A|i(_n(r — r'), one may then 
proceed as described in Sec. Ill and regard the projective 
variables f||(r, ui) and fj_(r, u>) as being two independent 
sets of bosonic variables. 

Let us briefly make contact with the quantization 
scheme described in Refs. [19, 21], where j (r,u;) is de- 
composed according to 



j N (r, w) = J Ne (r, w) + J Nm (r,w), 



with 



J-NeO^) 



he 



yw 2 Im£(w) f e (r, ui) , 



(89) 



(90) 



J Nl »( r > w ) = Ti y—J V x [\Z-Im«Mfm(r,w)]. 

(91) 

The connection between Eqs. (86)-(88) and Eqs. (89)- 
(91) is given by a gauge transformation (cf. Sec. Ill), 
which effectively redistributes the first term of Eq. (88). 
It is not difficult to prove that the total noise current 
as given by Eq. (89) satisfies the correct commutation 
relation (10) [with a(r, r',u>) from Eq. (77) together 
with Eqs. (82) and (83)] if f e (r, ui) and f m (r, w) are 
regarded as two independent sets of bosonic variables. 
Since J Ne (r,w) and J Nm (r,w) can be linearly related to 
j N ,,(r, ui) andj N± (r, w) and thus toj N (r, ui), the variables 

f e (r, ui) and f TO (r, ui) may be viewed as projective vari- 
ables, in the sense outlined at the end of Sec. III. Since 
Eq. (89) [with Eqs. (90) and (91)] is a separation of the 
noise current into a part attributed to a dielectric polar- 
ization and a part attributed to a (paramagnetic) mag- 
netization, the quantization scheme based on Eqs. (89)- 
(91), with f e (r, ui) and f m (r, ui) being bosonic variables, 
may be thought of as following from the general quantiza- 
tion scheme in the case where magneto-electric crossing 
effects can be a priori excluded from consideration. 



3. Local limit: other kinds of media 

The transition to the local limit is not a unique pro- 
cedure in general. Various kinds of locally responding 
(homogeneous) media, including non-isotropic ones, may 
therefore be derived as limiting cases from Eq. (57) . To il- 
lustrate this, let us represent a (k, ui) as given in Eq. (58) 
in a different orthonormal basis, where the new expansion 
will be non-diagonal in general, 

3 

<r(k,o;)= cTjj(k,w)e»(k,u;)e!-(k,u;). (92) 

The new basis vectors §j(k, ui) are related to the ones 
appearing in Eq. (58) by a unitary transformation, 



ej(k,w) = ^ ^fe(k,w)e fc (k,w), 



k=l 



U ik (\t, ui) = §i(k, lu) ■ e£(k, u). 



(93) 
(94) 



We may always choose the §j(k, ui) so that they are in- 
dependent of k, ej(k,w) i > ei(uj). If this choice can be 
made such that the new expansion coefficients, 

3 

aij(k,u)= c/ 4( k ^) cr fci(k,w)C/,7(k,w), (95) 

k,l=l 

may be approximately replaced under the k-integral ac- 
cording to 



<7jj(k,u;) i-> aij(\i.—>0,uj) =aij(uj) 



(96) 



when Eq. (92) is inserted in Eq. (57), then in this way 
the type of locally responding (homogeneous) anisotropic 
medium defined by Eq. (47) is recovered. [Equation (48) 
is then obtained by diagonalizing dij(w) by means of yet 
another (k-independent) unitary transformation.] Simi- 
larly, if the approximation (96) is generalized to include 
further terms of an (assumed) expansion of <7y(k, to) 
at k = 0, then quasi-local approximations of Eq. (57) 
are generated, by inserting the truncated expansion into 
Eq. (57) and integrating term by term to yield a linear 
combination of various derivatives of (^-functions. In pur- 
suing such approximation procedures — whose validity is 
to be examined in each case and which depends crucially 
on the choice of the transformation (93), (94) (i.e., of the 
new basis vectors) — it must be kept in mind that any 
approximate form of <r(r, r',u>) so derived has to con- 
form to all the general requirements on Q(r, r',u>). In 
fact, the most general kind of linear medium discussed in 
terms of local constitutive relations in the literature, the 
so-called bi-anisotropic medium (see, e.g., Ref. [23]), may 
be viewed in this way as a quasi-local approximation of 

<r(r, t',uj) [and thereby of Q(r, t',u>)] that incorporates 
derivatives of <5-functions up to the second order. 
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C. Spatially dispersive inhomogeneous media 

As already mentioned, knowledge of the medium prop- 
erties, i.e., of a (r, r', lu), is required in order to solve 
the eigenvalue problem (19) and perform explicitly the 
quantization the medium-assisted electromagnetic field — 
a task which, in general, cannot be accomplished in closed 
form. Nevertheless, to provide some analytical insight 
into the problem, let us consider media that combine the 
features of the media considered in Sees. IV A and IV B 
in an approximate fashion. 

1. Model 

We assume that the medium permits one to clearly dis- 
tinguish between the length scales associated with spatial 
dispersion and inhomogeneity, with the former scale be- 
ing sufficiently small as compared with the latter one. In 
this case, the medium can be regarded as having locally 
the properties of bulk material, and a(r,r',w) may be 
approximated as 

L i=l k 

x # L (r)# L (r> 4L kMe* Lk H, (97) 

from which the eigenfunctions of the associated operator 
are seen to be 

F lLk (r, w) - n-^W e 4k - r e lLk H. (98) 

Here, the medium is thought of as being divided into unit 
cells of volume which form a Bravais-type lattice, the 
cut-off function 8j_, (r) is unity if r is in the cell of lattice 
vector L and zero otherwise, ejLk(w) are, for chosen L, 
k, and lu, a triplet (i — 1, 2, 3) of orthogonal unit vectors, 
and the wave vector k runs over the reciprocal lattice. 
Note that, for each cell L, 

3 

<7 L kM = £ CT iLk(^)eiLk(^)e* Lk (a>) (99) 
»=i 

corresponds to the diagonal form in the (k, lu) domain of 
o(r,r' ,w) for bulk material [cf. Eq. (58)]. 

The main features of a (r, r', lu) as given in Eq. (97) can 
be summarized as follows, (i) tr(r, r',w) is zero when- 
ever r and r' are not in the same cell, so that fi 1 / 3 de- 
termines the length scale on which spatial dispersion is 
at most observed, (ii) The dependence on L of <7Lk(w) 
[Eq. (99)] for an inhomogeneous medium introduces an 
L-dcpendence into Eq. (97) which should be sufficiently 
weak, so that noticeable violations of the translational 
invariance of a (r, r', lu) may occur only on a length scale 
that is large compared with fi 1 / 3 . Needless to say that 
the main features do not essentially change if 9l(t) is 



replaced by another — but qualitatively similar — cut-off 
function. 

Let us denote by L(r) the particular lattice vector 
whose cell contains the point r, so that L(r) plays the role 
of a coarse-grained position variable. With the notations 
0L(r)(r') >-> 0[L(r),r'] and <7 L(r)k (u;) ^ <7 k [L(r), w], 
Eq. (97) together with Eq. (99) can be rewritten as 

a(r,r',c) = 0[L(r'),r]a[L(r'),r-r',u;], (100) 

with 

a [L(r'), r " *>] = ^ £ ^[L(r'), w]e ik -< M '> . (101) 

k 

Note that for arbitrary (continuous) values s, the func- 
tion 9(s, r) can be regarded as being symmetric. Using 
Eq. (98), wc find that Eq. (23) takes the form 

K<tJ iU ) = ^^l£fk[L(r'), W ]e*-(-'), 
" k 

(102) 

where 

3 

K k [L(r) M=J2 [L(r) , w] e, k [L(r) , c]e* k [L(r) , lu] 

i=l 

(103) 

{criL(r)kM i-> a ik \L(r),w], e iL(r)k (w) ^ e ik [L(r), lu]}. 

It can be shown that Eq. (100) [with Eq. (101)] in- 
deed contains (and, in a sense, interpolates) the two 
limiting cases studied in Sees. IV A and IV B. For the 
proof, we observe that in the case of negligible spatial 
dispersion, the cell size can be shrunk to zero, fl — > 0, 
so that the lattice vectors take on continuous values, 
L(r) — » r. As the lattice becomes finer and finer, the 
reciprocal lattice becomes more and more coarse, and, 
for r and r' unequal but in the same cell, all the points 
of the reciprocal lattice with k ^ give rise to rapidly 
oscillating terms in Eq. (101). In the limit — ► 0, 
these terms oscillate infinitely rapidly and average to zero 
(when applying the operator associated with Eq. (100) 
[with Eq. (101)] to any reasonable function), so that 
they may be set equal to zero. Taking also into ac- 
count that 0[L(r'),r]/fi — > 5(r — r') in this limit, we 
see that Eq. (100) [with Eq. (101)] indeed approaches 
Eq. (47) for vanishing spatial dispersion [note the corre- 
spondences (Tjk=o[L(r) = r, w]=<7i(r, uj) and e ik=0 [L(r) = 
r,w]=e i (r,w) ]. 

On the other hand, in the limiting case of an in- 
finitely extended homogeneous medium, there is no Li- 
dependence of the medium properties so that we are free 
to increase the cell size indefinitely, £1 — > oo. Conse- 
quently, wc may let 0[L(r'),r] — > 1 in Eq. (100) and 
ak[L(r),w] -> a(k,uj), fT 1 £ k - (27r)- 3 /d 3 fc in 
Eq. (101), which reveals that Eq. (100) [with Eq. (101)] 
approaches Eq. (57) as expected. 



11 



2. Magnetodielectric media 

To quantize the electromagnetic field in an inhomo- 
geneous magnetodielectric medium specified in terms of 
e(r, ui) and n(r, ui) — fi (r, ui), let us consider a medium 
that is both sufficiently weakly inhomogeneous and suffi- 
ciently weakly spatially dispersive, so that fl in Eq. (101) 
can be chosen on a scale intermediate between the scales 
of spatial dispersion and inhomogeneity. We may then 
approximately let L(r) be a continuous variable, L(r) — >r 
in Eq. (100), and yet, at the same time, approximately 
treat the k-sum in Eq. (101) as an integral, so that 
Eq. (100) [with Eq. (101)] approximates to [0[L(r'),r] -> 
0(r',r)} 

a(r,r',ui) = ^jf) J d 3 fc a(r', k, oj)^-^. (104) 

For a medium that is locally of the type described by 
Eq. (69), we may set 

<r(r, k, ui) = (T|| (r, k,ui) I - k x 7(1", k,ui)I x k, (105) 
where 

7(r,fc,w) = [a ± (r,k,ui) - a\\ (r, k, uj)]/k 2 > 0. (106) 

Assuming that in the k-integral in Eq. (104), a»(r,k,uj) 
and 7(1*, k,uj) may be approximated, respectively, by 
well-defined (and unique) long-wavelength limits a\\ (r, u) 
= limfc^o f || ( r 7 k, ui) and 7(1*, ui) = lim^o 7( r , k, ui), the 
cut-off function 0(r', r) has — due to the rapid oscillations 
of the exponential for large |r — r'| — no effect [with re- 
gard to an application of the operator associated with 
a(r, r',u) from Eq. (104)] and can be dropped, and we 
obtain, as a generalization of Eq. (77), 

a(v,v',ui) = cry (r', ui) I5(r - r') 

-Vx[7(r',w)/i(r-r')]xv'. (107) 

With the identifications 

<T||(r,w) = eoujlme(r, ui), (108) 
7(r, ui) = — K luiK(r, uu)/u (109) 

[cf. Eqs. (82) and (83)], Eqs. (78) and (79) generalize to 

Q(r, r', ui) = -ie uo [e(r', uo) - 1] Ts(r - r') 

- V x {«o[l - n(r', u)] Ts{r - r')} x V' (110) 

and 

j(r,w) = -ieoui [s{r,u>) - l]E(r,w) 

+ k V x {[l-K(r,w)]B(r,w)}+j N (r,w), (111) 
respectively. 



Unfortunately, Eq. (84) does not generalize to 

K'(r,r',uj) =aj /2 (r',w)7(5(r-r') 

T!7 1/2 (r»Vx/(5(r-r') (wrong!), (112) 

as could have been suspected. Indeed, straightforward 
calculation shows that, for spatially varying permittivity 
and permeability, the kernel (112) does not solve Eq. (18) 
[with cr(r, r', oj) as given in Eq. (107)], which implies that 
j N (r,w) cannot be related to the variables f(r,w) as in 
Eq. (85), with e(u>) and k(w) being simply replaced with 
their inhomogeneous counterparts e(r, ui) and /t(r, ui), re- 
spectively. In order to obtain an explicit expression for 

the kernel K(r, r', ui) required in Eq. (17), one has instead 
to return to Eq. (18) and solve it with cr(r, r',u>) from 
Eq. (107) — a problem that is, however, very difficult to 
solve in general. Although this does not at all limit the 
practical applicability of the theory [since all one typi- 
cally has to know about K(t,t' ,ui) is that it satisfies its 
defining equation (18) for the chosen conductivity (110)], 
it may be useful to have at hand, at least an approximate 
form for weak inhomogeneity, such as (App. D) 

K(r, r', w) = (r\ ,2 {ui) F(5(r-r')T«7 1/2 (^) V x JJ(r-r') 
+ ±[a ll (r,uj)-a ll {ui)}M (r, r» 
+ iVx{[ 7 (r, W )- 7 H]VxM„(r,r»}, (113) 
where 

M„(r,r»= ff| 1/2 H/i(r-r') 

± i7 -1 / 2 (w) V x mo(r, r', ui) I 

+ a~ 1/2 (ui)V x m (r,r»/x V', (114) 

with m (r,r',w) = -(47r|r - r'l)" 1 e - k-r' |/ Q (^) ^ a ^ = 
[7(a>)/<7|| (w)] 1 / 2 > 0. If the lack of exact knowledge of 

K(r, r',ui) really happens to be an obstacle in an appli- 
cation, one can alternatively resort to the approach on 
the basis of Eqs. (89)-(91), by simply replacing therein 
e(ui) and n{ui) by e(r, ui) and k(t,u)), respectively. 

V. GREEN TENSOR CONSTRUCTION FOR 
SPATIALLY DISPERSIVE BODIES 

Practical application of the quantization scheme re- 
quires the solution of the classical problem of the deter- 
mination of the Green tensor G(r, r',ui) for a given con- 
ductivity tensor Q(r, r',ui). Typically, one has to deal 
with systems of bodies each of which can be regarded as 
being homogeneous in its respective interior region. The 
physical surfaces of the bodies, including the boundary 
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surfaces between adjacent bodies, may be said to be the 
particular space regions where the material properties 
differ significantly from the intra-body (bulk-material) 
properties. Physical surfaces are therefore not mathe- 
matical ones but more or less fuzzy boundary layers. For 
not too small bodies, however, they usually contain only 
a small fraction of the overall material, so that they may 
often be approximately replaced with (sharp) mathemat- 
ical surfaces, with the idealization that the intra-body 
bulk-material properties hold immediately beyond them. 



satisfies the reciprocity condition. Corrections to the 
dielectric approximation in the form of surface currents 
or, equivalently, boundary conditions might be required 
for systems such as needle-shaped bodies, thin films or 
the like, where the dielectric approximation can be in- 
sufficient. However, as such corrections cannot be con- 
vincingly justified without detailed (model) assumptions 
about the specific nature of the physical surfaces, we do 
not consider them in the following. 



A. Dielectric approximation 

The application of the point of view just outlined to 
(systems of internally homogeneous) spatially dispersive 
bodies is commonly referred to as the dielectric approxi- 
mation, for which it is obviously required that the charac- 
teristic length scale of spatial dispersion is small in com- 
parison with the typical linear extensions of the bodies. 
Hence, making use of the same notation as employed in 
Eq. (100), we assume the conductivity tensor of a system 
of spatially dispersing bodies, in the dielectric approxi- 
mation, to be of the form 

Q(r, r', uj) = 0[L(r'), r] Q[L(r'),r - r', uj], (115) 

with the vector L labeling now the bodies in place of the 
lattice cells in Sec. IV CI; L(r') singles out the particu- 
lar body that contains the position r', and the quantities 

Q[L, r — r', uj] are the bulk-material conductivity tensors 
ascribed to the various bodies. (Regions outside all ac- 
tual bodies, if any, are formally viewed as bodies in this 

notation.) Note that Q(r, r',uj) as given by Eq. (115) 



B. Integral equations and surface impedance 
method 

The dielectric approximation renders it possible to for- 
mulate integral equations from which the Green ten- 
sor may be then derived. Here we outline the surface 
impedance method, which, in connection with its appli- 
cation to the calculation of Casimir forces, has recently 
given rise to controversial discussions (see, e.g., Refs. [24- 
27]) concerning the range of validity and the approxima- 
tions involved in this method. 

To begin with, let us consider spatially dispersive ma- 
terial described by a conductivity tensor Q(r, t',uj) [not 
yet approximated by a form like Eq. (115)], and let 
V denote some space region of interest. Further, let 

G aux (r, r', uj) be the Green tensor for an auxiliary prob- 
lem to be specified yet. Then, if j(r, uj) is an arbitrarily 
chosen current density in the frequency domain, given 
inside and/or outside V, and if E(r, uj) and B(r, w) = 
(«w) _1 VxE(r, uj), respectively, are the (classical) electric 
and induction fields associated with this current density, 
the identity 



/ d 3 rE(r,w).Vxr sll5 (r,r» + K(|W / d 3 rE(r,w)- G aux (r, r', uj) - (ie Q oj) 1 d 3 s Q(r, s, uj) ■ G aux (s, r', uj) 
Jv Jv L Jv 



-I 



d 3 r j (r, uj) ■ G aux (r, r', uj) + do(r) [E(r, uj) x e„(r)] • T aux (r, r', uj) 



dv 



+ ^ / da(r)[B(r,w) xe„(r)]-G 

a i ix (r,r',w) (116) 

JdV 



r 



(being a Green-type formula) holds, as can be proven cor- at r on the surface dV of V. Further, 
rect by partial integration and employing the reciprocity 

property of Q(r,s,w). Here, da(r) and e„(r) denote the 
absolute value and the unit vector of the surface clement 



j v ( r > w ) = (w^) 1 



V x V x E(r,w) 



E(r,uj)-ifi uj / d 3 r'Q(r,r',w) -E(r',w) 



(117) 
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and 

T aux (r, r', w) = (i/i w) _1 V x G aux (r, r', w). (118) 

Note that j y (r, uj) does not agree with j(r, uj), in general, 

as the r'-integral in Eq. (117) extends only over V. 

*— * 

Let us now assume that Q(r,r',uj) can be treated in 
the dielectric approximation according to Eq. (115), at 
least for the medium in V. In this case, j y (r, uj) agrees 

with j(r, uj) inside V, and Q(r,s,uj) in Eq. (116) can 
be replaced with the bulk-medium conductivity tensor 

<5(Ly,r — s, w) attributed to the medium in V (Ly is 

the vector labeling region V). Hence, if G aux (r, r', uj) is 
required, for r e V (at least), to obey the equation 

V x V x G aux (r, r', w) - G aux (r, r', w) 

-¥ow / d 3 sQ(Ly,r - s,w) • G aux (s,r',w) 
Jv 

= ?5{r-r'), (119) 

then Eq. (116) simplifies to 

(zMo^rw^MLvy) = ^r^v,^) 

+ / da(r) [E(r, uj) x e„(r)] • f aux (r, r', w) 

+ Mo" 1 / da ( r ) [B(r, w) x e„(r)] • G aux (r, r', uj), 

(120) 

where 

E^(r,^) = inou; f dVj(r>)-G aux (r',r, W ) 

= i/i u dVG aux (r,r',w) -j(r',w). (121) 

Since the characteristic length of spatial dispersion is as- 
sumed to be sufficiently small as compared with the linear 
extensions of V, we may extend, with little error, the s- 
integral in Eq. (119) to the whole space and, therefore, 

we may identify G aux (r, r', uj) in Eq. (120) with the corre- 
sponding bulk-medium Green tensor as given in App. E. 
Note that Eq. (120) may be viewed as a statement of 
Huyghens' principle (for r' inside V) and of the extinc- 
tion theorem (for r' outside V), see Refs. [23, 28], which 
have been well known suitable starting point for 

field calculations on the basis of integral equation meth- 
ods, see, in particular, Refs. [3-5] where the Wiener-Hopf 
technique has been used to construct solutions for a par- 
ticular functional form of Q(Ly , r — s, uj). 

The method of surface impedance (see, e.g., Ref. [29]) 
consists in the assumption that a linear relation between 
the tangential field components of E(r, u) and B (r, uj) 



exists on the surface dV, 

e„(r) x E(r,w) 
= Mo" 1 / d a (r')[F(r,r',^)xe„(r')]-[B(r',^)xe„(r')] 

JdV 

(122) 

(r on dV), where the tensor Z(r, r', lu) is the dyadic sur- 
face impedance. Equivalently, one may consider the in- 
verted form 

Ho ^(r.w) x e„(r) 
= f da(r')[E(r» x e„(r')] • Y(r', r, uj) ■ F„(r), 

JdV 

(123) 

where I „(r) = / — e„(r)e„(r) is a tangential projec- 
tor, and Y(r, r', to) may be referred to as the surface ad- 
mittance. Once the relation (122) [or (123)] has been 
adopted, one may convert Eq. (120) (with r' e V) into 
an integral equation for B(r, uj) [or E(r, uj)] inside V, 
which can then be solved in terms of E*- 11 ^ (r, uj) and the 
surface impedance, without specifying the medium prop- 
erties outside V . 

At first glance, the method has some advantageous fea- 
tures so that it has been enjoying a reputation in the 
literature. One of these features is that the necessity 
to know the medium properties everywhere outside V is 
replaced, so to speak, by the necessity to know the sur- 
face impedance for dV, which is a great reduction at first 
glance. Another one is that the application of continuity 
conditions may be sidestepped to some extent. Both of 
these points are to be qualified, however, and must be 
seen in the context of the following remarks. 

First, since the current density j(r, uj) outside V does 
not contribute to Eq. (121), the solution for E(r,cj) can 
be unique only if j(r, uj) is located completely inside V; 
solutions to the homogeneous integral equation can then 
be excluded from consideration. The surface impedance 
method applied to the volume V can thus yield, for given 

surface impedance, the Green tensor G(r, v 1 ,uj) only for 
both r and r' located in V (which may be sufficient for 
many applications, e.g., in the calculation of dispersion 
forces). Second, Eqs. (122) and (123) implicitly demand 
that the tangential components of the electric and the in- 
duction field uniquely exist on the surface dV (i.e, they 
should be continuous across the surface), or else ambi- 
guities were encountered. Quasi-local approximations of 
the conductivity tensor should hence be excluded if they 
contain magnetic-like singular terms [as the second one 
in Eq. (110)]. 

Provided that Eq. (122) [or (123)] holds, the method 
requires knowledge of Z(y,y',uj) [or Y"(r,r',w)], which 
plays the role of an external input. If Z(y,y',uj) [or 
y(r, r', uj)] is not known from the very beginning, one can 
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try to determine it a posteriori from the solution found, 
by appropriately specifying the effect of the medium on 
the electromagnetic outside the space region under con- 
sideration as well. We proceed with an example where 
this line can be pursued explicitly on the basis of the 
dielectric approximation. 



C. Example: Non-magnetic, planar systems 

Let V denote the slab-like region between two parallel 
planes z = and z = d, and let j(r, ui) be located com- 
pletely inside V. In order to apply the surface impedance 
method to the region V, we assume that the medium in V 
is non-magnetic and can be treated in the dielectric ap- 
proximation. Introducing Fourier transforms according 
to [r = (p, z)] 

E(r,w) = J d 2 qe tq - p E(z,q,w) (124) 

and taking into account the lateral translational invari- 
ance of the system, we may write Eq. (123) in the Fourier 
domain as 

fio 1 B(^,q,w) x e z 

= (2tt) 2 { [E(d, q, lu) x e z ] ■ Y(d, z, -q, to) 

- [E(0, q, lu) x e z ] ■ Y(0, z,-q,cu)}-T z (125) 

[ I z = I — e z e z ] , where z takes on the two values z = and 
z — d (the two terms enter with opposite signs because 
of the opposite surface normals on the two sides of dV). 
Similarly, the Fourier transformed version of Eq. (120) 
reads 

E(z',q, W )0(Ly,z')-E (in) (^,q,^) 
= iLi Lu(2ir) 2 {[E(d, q,w) x e z ] ■ T aux (d, z 1 , -q,w) 

- [E(0, q, lu) x e z ] ■ r aux (0, z', -q, w) 
+ Mo 1 ©( d ' c L> w ) x e z] • G aux (d, z', -q,w) 

- / i o 1 [B(0,q, W )xeJ.L(0, Z '-q^)}. (126) 

Note that G aux (r, r', w) and T aux (r, r', w) are transla- 
tionally invariant, because they refer to bulk material. 
By means of Eq. (125), Eq. (126) takes the form 

E(z',q,c)0(L y ,z')-E (in) (^,q,^) 

= E(d, q, lu) ■ I z ■ R(d, z' , q, lu) 

-E(0,q,o;)- J, • fl(0,«',q, w), (127) 



with (C = 0, C = d) 

^(C,^',q,w) =^o^(27r) 2 e z x { 2', -q, w) 

+ (27r) 2 f (C, d, -q, w) • / , ■ Gi (d, z', -q, w) 
-(2^) 2 y > (C,0,-q,^). / z .Gi(0,« , ,-q,w)} (128) 

[note that e z ■ R((, z', q, lu) = 0]. 

From Eq. (127) it follows that (Appendix F) 

E(0,q,w)-7, = [E( in )(0+,q,w)-C» 

-E (in) (d-,q,w) .£)«] • - (129) 

E(d,q, W )./ 2 = [E( in )(0+,q,a;)-I» 

-E (in) (d-,q,w) • -5 •]?»]», (130) 

where 

2=2(q,w)= / z -[J+^(0,0+,q,w)]-/„ (131) 

B = B(q,u) = / z -^(0,d-,q,w)- J z , (132) 

£=C(q,w) = -J z -.R(d,0+,q,w)- J z , (133) 

D = .D(q,w) = J z -[J-fl(d,d-,q,w)]- F z , (134) 

and the superscript jj denotes matrix inversion with re- 
spect to the 7 z -space: A$ = I z ■ (I z ■ A ■ I z ) _1 • I z . 
The solution to Eq. (127) can now be readily obtained 
by substituting Eqs. (129) and (130) back in the right- 
hand side of Eq. (127). Specifying in Eq. (121) the 
source of E*- ln ' (r, lu) so as to correspond to a point dipole 
situated in V, j(r, lu) = (inou)~ 1 p5(r — r p ), from the 
corresponding position-space solution E(r,w)^(Ly,r), 
one can then read off, according to E(r, lu)6(Lv, r ) = 

p • G(r p , r, cj)8(Lv, r)6(Lv, r p ), the (interesting part of 
the) Green tensor expressed in terms of the surface ad- 

mittance, i.e., in terms of the quantities Y(z, z' , q, lu) (z 
= 0,d; z' = 0,d). 

In order to illustrate the calculation of the admittance, 
let us assume that the (non-magnetic) media to the left 
and right of z — and z = d, respectively, are (in the 
sense of the dielectric approximation) homogeneous half- 
spaces. Using the indices 0, 1, and 2 to distinguish the 
left (z < 0) half-space, the volume V (0 < z < d), and 
the right half-space (z>d), respectively, and associating 
with each of the three regions the corresponding bulk- 
medium Green tensor Gj(r, r',Lu) (App. E) and the as- 
sociated auxiliary tensor Tj(r, r', lu) defined according to 
Eq. (118), one can derive the admittance by applying the 
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continuity conditions. Introducing the abbreviations 

M(q,w) = -/ z -Gi(0,d-,q,w)- / z , (135) 

N(q, lo) = -T z - [Gi(0, 0+, q, w) + G o (0, 0-, q, w)] • 7 Z , 

(136) 



P(q,w) = / a -[Gi(d,d-,q,w) + G 2 (0,0+,q,w)]- /„ 

(137) 



S(q,w) = J* • Gi(d,0+,q,w) • J«, 
T(q,w) = F z • ri(0,d-,q,w) • J Z) 



(138) 
(139) 



tf(q,w) = /* • [ri(0,0+,q,w) + r o (0,0-,q,w)] • J z , 

(140) 

f(q,u;) = -7, • [ri(d,d-,q,w) + f 2 (0, 0+, q, u)] • /„ 

(141) 



W(q,w) = -J z -ri(d,0+,q,w)- I z , 



(142) 



one eventually finds (App. G) the Y(z, z', q, ui) (z = 0, d; 
z' = 0, d) to be (for notational convenience, the arguments 
q, lo are suppressed) 

^( '°) = 7^-v?P ■ ^ - P • P"] ■ [M ■ P* — iV ■ Sty, 



(2tt)2 



(143) 



Y(0,d) = -^[U-N* - T • M"] • [P • - S • A»] 



(144) 



1 



F(rf,0) = — — [F.p« - VF- 5 J ] • [M • P" — TV ■ S 1 *]*, 
( 27r ) 

(145) 

Y(d,d) = [V^- M tt - W • A tt ] • [P • i? s - 5 • ivty, 

(146) 

which are, quite naturally, far more complicated than the 
simple Kliever-Fuchs expressions mentioned in App. E. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have developed a rather general quantization 
scheme for the macroscopic electromagnetic field in arbi- 
trary linearly responding media, which offers a unified ap- 
proach to QED in linear media. Describing the medium 
response by a non-local conductivity tensor, any of the 
possible electromagnetic features of a linear medium in 
equilibrium is covered by the scheme, in particular, spa- 
tial dispersion. Central quantities of the scheme are the 
noise current that is intimately connected with the ab- 
sorption necessarily observed in any linear medium in 



equilibrium, the bosonic dynamical variables associated 
with the noise current, and the Green tensor of the phe- 
nomenological Maxwell equations, in which the medium 
properties enter via the conductivity tensor. 

From a careful analysis of the dynamical variables and 
(quasi-)local limiting forms of the non-local conductivity 
tensor, we have shown how quantization schemes pre- 
viously developed for locally responding media can be 
recovered as special applications of the general quantiza- 
tion scheme. In particular, a locally responding magne- 
todiclectric medium can be viewed as a special quasi-local 
limiting case of an isotropic, spatially dispersive medium 
without optical activity, where the (local) dielectric per- 
mittivity and magnetic permeability arc just two con- 
tributions to one and the same quasi-local conductivity 
tensor. As a result, application of the general quantiza- 
tion scheme shows that the electromagnetic field in such a 
medium can be quantized by using a single set of bosonic 
variables. 

Generally, the use of a single set of bosonic variables 
means that the noise current which enters the macro- 
scopic Maxwell equations is not divided into parts (as- 
sociated, e.g., with a polarization and a magnetization) 
regarded as representing independent degrees of freedom, 
but is rather treated as an entity. This may be particu- 
larly advantageous for future studies of (quantum) elec- 
trodynamics in moving media, simplifying the discussion 
of transformations to different frames of reference. How- 
ever, the theory also admits, by appropriate projection, 
the use of several independent sets of bosonic variables, 
which in fact corresponds to the neglect of certain kinds 
of interactions in the sense of super-selection rules. 

Since exact solutions of Maxwell's equations are not 
available in closed form in general, even more so if spa- 
tial dispersion is taken into account, one has to resort 
to approximation methods to obtain explicit expressions 
for the Green tensor, the latter being one of the corner- 
stones of the theory. To assist in such intentions we have 
considered in some detail the dielectric approximation 
for the conductivity tensor, which consists in approxi- 
mating the conductivity tensor of a system of spatially 
dispersive bodies by joining together bulk-medium con- 
ductivity tensors (which are routinely handled in recip- 
rocal space) . Although the information relevant to phys- 
ical surface regions is lost in this way, the dielectric ap- 
proximation has been a key tool to render tractable elec- 
tromagnetic propagation problems in spatially dispersive 
media. Accepting the dielectric approximation, the prob- 
lem of finding the Green tensor becomes then solvable via 
integral-equation and surface-impedance techniques. 

As already mentioned, diamagnctic media are not 
covered by the quantization scheme developed in this 
paper — a scheme that exhausts the possibilities offered 
by the linear-response framework. Furthermore, it should 
be pointed out that the scheme does also not automati- 
cally apply to linearly amplifying media. Although both 
types of media do not really fit into the linear-response 
framework, they may be forced into it, but not without 
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reservations and alterations of the whole scheme. Clearly, 
the concept of linear amplification has a range of applica- 
bility very much smaller than that of linear dissipation. 

Concluding, this work provides the most general quan- 
tization scheme for the electromagnetic field in linearly 
responding, absorbing materials to date, from which pre- 
viously given schemes can be recovered as limiting cases. 
It will serve as a foundation for investigations of surface 
plasmon effects involving strong spatial dispersion, and 
as a starting point for the investigation of moving media. 

APPENDIX A: DERIVATION OF EQ. (8) 



From Eq. (A2) it is seen that 
ImH(r, r',u) = 

T 2 

--^ST- /£(r-r')-MoRe[wQ(r,r>)]. (A8) 

Insertion of Eq. (A8) in Eq. (A7) and restriction to real 
frequencies leads, upon recalling Eq. (2), to Eq. (8). 

APPENDIX B: NON-UNIQUENESS OF THE 
KERNEL K(r,r',u) 



The linear integro-differential equation (8) can be rep- 
resented as 



/ 



d 3 stf(r,s,w) • G(s,r',w) = 16{r-r'), (Al) 



where the integral kernel 

H(r, r', oj) = V x V x ?S(r - r') 



The transition from K(r,r',oj) to K' (r,r' ,u>) accord- 
ing to Eq. (24) can be re-interpreted as a redefinition of 
the dynamical variables f(r,w) and V(r,to) according to 

f (r, u) = J dV V(r, r', u) ■ f'(r', w ), (Bl) 

f t (r,w)= fdVf^r'.wj.f'^r'.w). (B2) 



-^I5(t- r') - i/i wQ(r,r» (A2) Inserting Eq. (Bl) into Eq. (17) yields 



is reciprocal, 



H(r,r',uj) = H T (r',r,uj), 



(A3) 



since Q(r,r',w) is reciprocal. Hence, the transposed 
equation of Eq. (Al) takes the form 

J d 3 s G T (s, r, w) • H(s, r', u>) = ?S(r - r'). (A4) 

Multiplying from the right with G(r',s',w), integrating 
over r', and using Eq. (Al), one can see that the Green 
tensor is also reciprocal, 



G(r,r» = G T (r',r,c). 



(A5) 



Because of Eq. (A5), the complex conjugate of Eq. (A4) 
reads 



J d 3 sG*( 



r,s,w) -H*(s,r',uj) = IS(t-t'). (A6) 

Taking the dot product of Eq. (Al) from the left with 

G* (s', r, u>) and integrating over r, taking the dot product 

of Eq. (A6) from the right with G(r',s') and integrating 
over r', and subtracting the two resulting equations, one 
derives 

ImG(r, r',u) = 

I / d3s '^ r ' s ' w )'[ Im ^( s ' s '' w )]'^*( s '' r '' w )- 

(A7) 



j x (r.^) = / d :i r'A"(r.r'..-).f'(y..-). (D3.) 



where K'(r, r',w) is just given by Eq. (24). With re- 
gard to the transformation (Bl) and (B2), the signif- 
icance of replacing Eq. (25) with Eq. (26) is that the 
variables f'(r, w) and f't(r,u;) are uniquely expressible 
in terms of the f(r,w) and F(r, w), and so are on an 
equal footing with them — the unitary operator associated 

with V(r, r', ui) uniquely maps a set of bosonic variables 
onto a fully equivalent set of bosonic variables. Hence, 

V(r, r',oj) may be thought of as being included in the 
chosen set of dynamical variables. In this sense, it is suffi- 
cient to base the calculations in Sec. Ill on the Hermitian 
operator associated with the integral kernel K(r, r', ui) as 
defined by Eq. (23). 

It is worth noting that the operator associated with 

K'(r,r,u)) as defined by Eq. (24) is non-Hermitian 

whenever V(r, r',w) is non-trivial. To see this, let us 
conversely assume that the operator associated with 

K'(r, r, co) is Hermitian, 

J d 3 sK(r,s,u)-V{s,r',oj) 

= j d 3 sV+(s,r,uj) ■ K(s,r',uj). (B4) 

Applying from the left the operator associated with V 
and from the right the operator associated with V + and 
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recalling Eq. (25), one sees that 



J d 3 sV{r,s,u))-K{s,r',Lu) 

= J d 3 sK(r,s,uj) • V*+(r',s,w). (B5) 



Applying the operator associated with K from the left 
to Eq. (B4) and from the right to Eq. (B5) and com- 
paring the results, one finds that the operators associ- 
ated with V and a commute [recall Eq. (18)], so that 
the operator associated with V maps each (possibly de- 
generate) eigenspace of the operator associated with a 
onto itself. Specifically, this implies that the operator as- 
sociated with V commutes with the spectral projectors 
of the operator associated with a [and, therefore, also 
with the projectors (34)]. Since the spectral projectors 

of the operators associated with a and K are the same 
[cf. Eqs. (22) and (23)], the operators associated with V 
and K also commute. But then, since the operator as- 
sociated with K is invertible (being a positive operator), 
Eq. (B4) [or Eq. (B5)] shows that the operator associated 

with V is Hcrmitian, i.e., 



V(r,r',a,) = V+(r',r, W ). 



(B6) 



Since the operator associated with V is also unitary, in 
the diagonal expansion 



J dav(a, 



w)F(a,r,w)F*(a,r , ,w), (B7) 



one must have v(a,u>) = ±1 for each a, which means 

that the Hermitian operator associated with K'(r, r', ui) 

can differ from the operator associated with K(r,r',w) 
only by the trivial type of unitary transformation that 
merely replaces some of the basis functions F(a, r, ui) 
with — F(a,r, ui). Conversely, this shows that any (in 

this sense) non-trivial V(r, r', u>) necessarily yields a non- 

Hcrmitian K'(r, r', w). 



APPENDIX C: REDUCED STATE SPACE AND 
SUPER-SELECTION RULE 



Let us consider the state space spanned by the Fock 
states associated with f(r, ui) and V(r,u>) so that an ar- 
bitrary, normalizable state \<p) in this space can be rep- 



resented in the form 
|0) = |0> <O|0> 



+ y~] / dwi / d 3 ri fcl (ri,^i) |l fcl (ri,wi)) 
kl =i J " J 

+ V / dwi / dcj 2 / d 3 ri / d 3 r 2 

> u , J0 JO J J 

x faik* (ri,wi,r 2 ,w 2 ) |l fcl (ri,wi),l fe2 (r2,w 2 )) 

+ ..., (CI) 



where 

/ fc (r,w)|0>=0, (C2) 
/>, W )|0) = |l fc (r,o;)>, (C3) 
fl N (r N ,w N )---fl i (r 1 ,w 1 )\0) 

= |l fcl (ri,wi), . . . , l kN (r N , uj n )) . (C4) 

The normalization of \<fi) can be obtained by using the 
formula (which can be viewed as a special case of the 
Bloch-De Dominicis theorem [14]) 

(0| fk M (rif, wm) • • • fk! (ri, wi) 



x 



ft (v' 1 ,u;[)...fl,(v' N ,uj' N )\0) 



( (0 1 0) = 1; Sn, group of permutations of N objects). 

In order to construct a reduced state space in which 
the operators f\(r, ui) and fj(r',a/) defined by Eq. (35) 
behave like bosonic operators, let us first introduce states 
0) A according to 

Ai(r,w)|0) A = 0, (C6) 
f Xi (T, ui) |0) A , = |0) v / Ai (r, W ) (A A') (C7) 

(a(0|0)a = 1), such that 

A 

'"Va- 



|0)=(g) 



(C8) 



A=l 



Now let us introduce, for each A, an orthogonal projector 
P\ as the sum of orthogonal projectors P^ N \ 



(N) 
X ' 



N=0 



pWt = p(w) 



and specify in such a way that, when applied to 

a quantum state of the form (CI), it picks out the 



(C9) 

(CIO) 
(Cll) 
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(N + l)th term on the right-hand side of Eq. (CI) and in- 
corporates N position-space projection kernels belonging 
to the chosen value of A, 

- |0> A (0| A (C12) 



i pOO p poo p poo p 

^I^l 



X /It 1 ( r l. w l)/{fe( r 2. u 2) ' ••/Afc JV ( r JV,^v)-F , A VAfcjvlriV^iv)^-!^-!,^-!)---/^^^!,^!) (C13) 

I 



(N = 1, 2, . . .)• It is not difficult to prove that Eqs. (CIO) 
and (Cll) are fulfilled, where the latter equation fixes 
the normalization factor 1/N\ in Eq. (C13), and that, in 
view of Eqs. (C8) and (Cll), the commutation relation 



(C14) 



holds. 

We may now define a reduced state space that contains 
only those (normalizable) vectors that have the separable 
form 



l^> (rcd) -(g)l^ 

A=l 

A \<t>) X = \<t>) A , 



A ' 



(C15) 
(C16) 



with each vector \<fi) x being, by construction, a superpo- 
sition of vectors 



\N) X = J2 dw 1 d 3 ri ...J2 dw N /d ; 

fc! J ° J k N J ° J 



r N 



x CxkL-.XkN (rij wi, • ■ ■ , rjv, wn) 
x |l Afcl (n,u;i), . . . , lAfe,v(rjv,wjv)) , (C17) 

where, in analogy to Eq. (C4), 

lUfei (ri,u;i), . . . , lAfc N ( r iv,wjv)) 

= /l feiv (riv,Wiv)---/l fel (ri,wi)|0) A . (C18) 

The important feature of these states is that the result of 
performing the integrations in Eq. (C17) is not changed 
if the wave function C xkl ... xkN (ri, u>i, ■ ■ ■ , rjv, lon) is re- 
placed according to 

Cxk!...Xk N (^i, wi, • • • , rjv, wat) 

^ / d3r '"' / ^ ^) fc i fc i( r i' r i' w i)-" 

x {P\)k N k' N {rN,r' N , uJN)C X k[...\k' N (r'i, wi, • • •, r^, Wat). 

(C19) 



It is also not changed if CAfci...Afc N (ri, Wi, ■ ■ ■ ,r N ,w N ) is 
symmetrized with repect to the labels 1, . . . , N. Wave 
functions that can be reduced to the same standardized 
wave function by these operations are thus fully equiva- 
lent representatives of the same vector. Without loss of 
generality, one can thus adopt the convention to employ 
only such standardized wave functions. 

The commutation relation (36) implies that 



' ,,/^- £ ^fc(r,w) 



= f{, k , (r', J) + e5 xx , {P x ) kk , (r, r', u)6(uj - J) (C20) 

with e being a parameter. As Eq. (C20) is a similarity 
transformation, it generalizes to 

e e/Afc(r ^ ) F[/], fe/ (r',w')]e _£/Afc(r '" ) 

= F[fl k , (r', J) + eS xx , (P x ) kk , (r, r» J(w - w')} 

(C21) 



of /j,,, (r',u/). Comparison of the terms of first order in 



where F = F[f x , k , (r', is any well-behaved functional 
e on both sides yields 

~f xk (r,u),F[fl k ,(r',uj')} 



F[flk'(*'^') 



+ e6 xx ,(P x )k k '(r,r',u ) )5(u>-u>')} \ . (C22) 

J £ = 

Let us consider the particular functional F]^[fy k ,(r', u) 1 )] 
appearing in Eqs. (C17), (C18), 

F N [fl, k ,(v',J)] 

POO P POO p 

= V / dwi / d 3 n . . . V / duj N / d 3 rjv 
fel Jo J k N Jo J 

x C xkl ... X k N (ri,uJi, ■ ■ ■ ,r N ,uj N ) 

></l fe >A^;v)--71 fel (ri,wi). (C23) 
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If the convention to use only standardized wave functions 
is adopted, one may write 

F N [fl k , (r', J) + e5\\> (P x ) kk , (r, r', J)5{u - J)] 

(C24) 

which means that the right-hand side of Eq. (C22) may 
be evaluated, for this functional, just as an ordinary func- 
tional derivative, i.e., 



f xk (v,cj),F N [fl >k> (r',u;') 



5F N [fl k ,(r',u')} 



(C25) 



But since, due to the definition of the reduced state space, 
only commutators of the type (C25) (for all N) are re- 
quired, and since Eq. (C25) can be obtained from Eq. (43) 
in the same way that Eq. (C22) has been obtained from 
Eq. (36), Eq. (43) is generally valid for the reduced state 
space. 



APPENDIX D: DERIVATION OF EQ. (113) 

For notational convenience, let us write here the inte- 
gral equation (18) in the compact operator form 



KKJ = a, 



with /C, and a being, respectively, the operators associ- 
ated (for chosen uS) with the integral kernels K(r,r',co) 
and <r(r, r',w). Accordingly, Eqs. (24) and (26) read 
K! =JCV and V^V — VV^ =T, respectively (1, unit op- 
erator). Assuming that the (Hcrmitian and positive) op- 
erator a takes the form 



a = cr + eoi, 



(D2) 



with e being a small, real parameter, we may try to find 
a solution to Eq. (Dl) by the perturbative ansatz 



JC = IC + elC[ + . . . , 



(D3) 



where /C is a solution to Eq. (Dl) for e = 0. Substituting 
Eqs. (D2) and (D3) into Eq. (Dl), we sec that the first- 
order correction tC[ obeys the equation 



ICqIC'i + K-iICq = <ti, 



(D4) 



which determines the Hermitian part of K,'olC'\ (recall 
that o i is Hcrmitian), whereas the anti-Hermitian part 
is left undetermined. The solution to Eq. (D4) may there- 
fore be written as 



£i = \{a 1 +A)M , 



(D5) 



where AA$ = (/C ) 1 , and A = — A* is an arbitrary anti- 
Hermitian operator, which may be simply set to zero. 



Note that the freedom to choose A corresponds to the 
freedom to choose V. We hence obtain to first order in 
a- cr 



/C' = /C + i(a-a )Xo. 



(D6) 



If do, and a are identified with the operators associ- 
ated with the kernels (77), (84) and (107), respectively, 
Eq. (D6) is just the operator equivalent of Eq. (113) 

[K'(v,v',uj)~K(ry,uj)}. 

To calculate explicitly the integral kernel M (r, r',iv) 
associated with Mo, we consider the Fourier representa- 
tion 



£<5(r,r» - j^f / d3fc ^o(k,^> 



) e ik.(r-r') 5 ( D7 ) 

where, according to Eq. (71), 

%w)=(jj /2 (w)[/iQ(a))kxF] (D8) 

[a(uj) — [y^/o^oj)] 1 / 2 > 0], which shows that the kernel 
M (r, r', u) can be given by the Fourier integral 

/ \ -1/2/ 



M (r,r',Lu)=a- 1/z (co 



, (D9) 



(Dl) an d h is n °t difficult to prove that 



{oj) 2 k 2 

, , o 1^ x 1 x k . . _ , 

+ «M 2 7—7^12- (D10) 



l + a{uj) 2 k 2 
Introducing the function 

1 f gik-(r-r') 

m (r,r ,u) = ——- / d d k 

(ziry J a(Lu) 1 + k z 

= -(47r|r - r '|)-i e -l r - r 'l/«(^), (Dll) 



we may rewrite Eq. (D9) [with Eq. (D10)] to ob- 
tain Eq. (114). Note that the (Yukawa-type) function 
TO (r,r',cj) satisfies the equation 

[-A + a{u)- 2 ] m (r, r', w) - S(r - r') (D12) 

together with the boundary condition m (r, r',u>) — > 
for |r — r'l — > oo. 



APPENDIX E: BULK-MEDIUM GREEN 
TENSOR AND KLIEVER-FUCHS IMPEDANCE 

For (translationally invariant) bulk material, 
d 3 fc 



,r>) = /(2^ elk ' ( " r,) ^ (k ' w) ' (E1) 



20 



the solution to Eq. (8) has the form 



G (0) (r,r» 



d 3 fc 
(2tt) 3 



eik .(r-r') G (0) (kjW); (E2) 



where (k, uj) is the solution to a simple 3x3 matrix 
equation. In particular, for an isotropic medium without 
optical activity, 



< — > kk kk 

Q(k,w) = Qn(fc,w)p- + Q ± (fc,w) ( / - 



fc 2 



one finds that 



G (0) (k,w) 



7 - kk/fc 2 kk/k 2 



D±(k,to) Du(k,U>) 



where 



D_i_(k, lu) = k 2 — lu 2 /c 2 — ifioujQ_i_(k, uj), 
D\\(k,w) = lu 2 /c 2 + ifj, u>Q\\(k,uj), 



(E3) 



(E4) 



(E5) 
(E6) 



and 



Q\\(±)(k,uj) = -i£ov[e\\( ±) (k,u}) - 1] (E7) 
in 'dielectric' notation. In the general case 
kk 



Q(k,w) = Q||(k,w) fc2 

kk\ ~ , x /<- kk 
•0±(k^)-M-TJ 



fc 2 



fc 2 



(E8) 



Z>H (fe, cf) changes to D||(k, w), D±(k,uj) changes to the 
tensor 

D ± (k,uj) = k 2 - ui 2 /c 2 - i/j, ujQ ± (k,uj), (E9) 

and the first term on the right-hand side of Eq. (E4) has 
to be replaced according to (7fc = 7 — kk/k 2 ) 



I - kk/k 2 

D±(k,u) 



I k ■ D±(k,uj) ■ I k 



Ik- (E10) 



It may be convenient — particularly with regard to sys- 
tems that are translationally invariant only in a plane, 
say the xy plane — to rewrite Eq. (E2) as [r = (p, z), k = 
(q,/3)] 

G(°)(r,r» = J d 2 qe lcl < p ~ p ">G^{z,z , ,q,Lu), (Ell) 
where 



G {a \k,u). (E12) 



The analytical properties of the integrand in Eq. (E12) 
with respect to (3 depend on the specific decay to 



zero of Q(r — r',w) for |r — r' — > oo. For sufficiently 
rapid decay Eq. (E12) admits an evaluation by con- 
tour integration in the complex (3 plane, which will be 
governed — focusing again on isotropic media without op- 
tical activity — by the solutions f3 = (q, a) of the two 
dispersion equations D±^(k,oj) = 0. In contrast to the 
case of spatially non-dispersive material, these equations 
are transcendental with respect to (3 rather than polyno- 
mial, so that nothing can be said about the number of 
their solutions in general. Specifically, if there are more 
than two functions f3„(q,u>) [and/or one or more func- 
tions /?|l(<7, u>)], the medium is said to support 'additional' 
(inhomogeneous plane-)waves. 

We close this appendix with the following (perhaps not 
well-known) observation. Inserting Eq. (E4) in Eq. (E12) 
and setting z — z' = 0±, and making use of the decompo- 
sitions 

7 - kk/k 2 = e s (q)e s (q) +e p (k)e p (k) (E13) 
7 Z = e s (q)e s (q)+qq/< Z 2 , (E14) 

with e s (q) = qx e z /q and e p (k) = k x q x e z /kq being 
polarization unit vectors, one can show that 

F z -G^(0±,0,q,^)- 7 Z 

= (i^y 1 [Z s (q,uj)e s (q L )e s (q) + Z p (q,uj) qq/q 2 ] , 

(E15) 



where 



1/IqUJ 



d(3 



o */30± 



(2tt)3 D ± (k,u) ' 



Z p (q,u) = inoui 

, (2tt) 3 k 2 



D±(k,uj) D»(k,u)) 



(E16) 



• (E17) 



With Eqs. (E5)-(E7), Eqs. (E16) and (E17) are recog- 
nized (up to a trivial factor) as the surface impedance 
expressions first derived by Kliever and Fuchs [30] for a 
spatially dispersive half-space by assuming specular elec- 
tron reflection (see also Ref. [27]). 



APPENDIX F: DERIVATION OF EQS. (129) AND 

(130) 

From inspection of Eq. (127) it is seen that only the 
tangential components of the electric field, which are as- 
sumed to be continuous at z' = and z' = d, contribute 
to the right-hand side of this equation. Wc therefore 
evaluate Eq. (127) at z' = 0+ and z' = d— and take the 

tangential (7 z ) component thereof to obtain the linear 
equations (the arguments q and u> are kept fixed in this 
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appendix and are suppressed for notational convenience) 

E(0+)- /*■[/+ £(0,0+)]- F z 
- E(d-) • • £(d, 0+) • F z = (0+) • I z (Fl) 

and 

E(0+) • Y z ■ R(0, d-) ■ Y z 

+ E(d-) • F z ■ [F- R(d,d-)] ■ 7 Z = E (in) (d-) • F z , 

(F2) 

respectively, which are to be solved for the tangential 

components E(0+) • I z and E(d— ) • 7 Z . To represent 
the solution in a compact form, we assign to any matrix 

A satisfying I z ■ A ■ I z = A its inverse on the / z -space, 

A$= I z ■ (I z ■ A ■ I z)^ 1 ■ I z . It is straightforward to see 
that the block-matrix formula 



together with Eqs. (131)-(134). Recalling again the con- 
tinuity of the tangential component of the electric field, 
we are left with Eq. (129) and (130). It should be 
noted that the above inversion procedure fails at partic- 
ular values of the (suppressed) arguments q and cj, be- 
cause of singularities. However, as we are dealing with a 
lossy system, such singularities — corresponding to guided 
waves — may appear only when lmw<0 (for real values 
Ofq). 



APPENDIX G: DERIVATION OF 
EQS. (143)-(146) 



A B\ (Ct-[A-C*-B-Dt]t A*-[C-A*-D-B*} 
K C DJ \pt-[B-B-A-C*]l B-[D-B-C-A$] 

V o i z 



(F3) 



is generally valid whenever the requisite inverse elements 
exist, so that the solution to Eqs. (Fl) and (F2) can be 
written in the form 

E(0+)-F z = [E (in) (0+)-C tt 

-E^id-) -B] ■ [A-C* - B-D*]*, 

(F4) 

E(d-)- T z = [E (in) (0+)-I» 

-E (in) (d-) -B] ■ [C-A* -D-B]* 

(F5) 



Let us attribute to the three regions j = (z < 0), 
j = 1 (0 < z < d), j — 2 (z > d) bulk-medium con- 
ductivities Qj{v — r',w), which combine to the over- 
all conductivity tensor in the sense of Eq. (115) [j <-> 
L]. For each of these regions, we construct, according 
to Eq. (118), the translationally invariant bulk-medium 

Green tensor Gj(r, r',u>) and the associated auxiliary 

tensor T j(r,r' ,u>) in terms of their Fourier components 

Gj(z,z',q,u)) and F j(z, z' ,q,cu) [defined according to 
Eqs. (Ell), (E12)]. The three regions are thus described 
on an equal footing so that for the field in each region, 
an equation similar to Eq. (126) holds. Evaluating the 
tangential components of these three equations (which 
together determine the field in all of space) and using 
the continuity conditions at z — and z = d, one obtains 
two sets of equations for the tangential boundary values, 



{[E(d) x e x ] • Ti(d,d-) + [Mo 'B(d) x e z ] • Gi(d.d-) - [E(0) x e z ] ■ ri(0,d-) - [^o ^(0) x e z ] • Gi(0,d-)} • I x 

= -{[E(d) xe 2 ]-f 2 (0,0+) + [Mo 1 B(d) x e z ] • G 2 (0, 0+)} • T z (Gl) 

and 

{[E(d) x e z ] ■ ri(d,0+) + K'Bfrf) x e z ] ■ Gi(d,0+) - [E(0) x e z ] • ft (0,0+) - [^^(O) x e z ] ■ <ft(0,0+)} • 7 Z 

= {[E(0) x e z ] ■ F o (0, 0-) + [mo '6(0) x e z ] ■ G (0,0-)} • T z , (G2) 
I 

where, for notational convenience, the arguments^ and f\ have been suppressed. Solving these linear relations 
u of E and B, and the arguments -q and uj of Gj and for B(0) x e z and B(d) x e z in terms of E(0) x e z and 



22 



E(d) x e z , we need only compare the result with Eq. (125) 
to verify Eqs. (143)-(146). This can be conveniently done 
by representing Eqs. (Gl) and (G2) in the form 



obtains 



„-iB(0)xe 
M I B(d) x e 




E(0) x e, 
E(d) x e, 



- - ■ 

\V WJ 

(G3) 

where M, N, P, S, T, U, V, W are given in Eqs. (135)- 
(142). Rewriting Eq. (125) (z — 0, d) in an analogous form 
and applying to Eq. (G3) the inversion formula (F3), one 




-J (0,0) -Y(0,d) \ = J_ 
Y(d,0) Y(d,d) J ( 2 < 



< IPt-[M-Pt-N-St]* Mt-[P-Mt-S-Nt]t\ 
K \j}*-[N-S»-M-pl]i Nt-[S-&-P-Mt]*)' 



which immediately leads to Eqs. (143)-(146). 
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